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Abstract. Galaxy clusters are hotter and more X-ray luminous than in quiescence while they undergo major mergers, which 
also transiently increase their strong-lensing efficiency. We use semi-analytic models for both effects to study how cluster 
dynamics in different dark-energy models affects the X-ray selected cluster population and its strong-lensing optical depth. We 
find that mergers increase the number of observable X-ray clusters by factors of a few and considerably broaden their redshift 
distribution. Strong-lensing optical depths are increased by a very similar amount. Quite independent of cosmology, X-ray 
bright clusters above a flux limit of 10~ 13 - 5 ergs -1 cirT 2 produce ~ 60% of the strong-lensing optical depth, and only ~ 1% 
above a flux limit of 10" 11 5 erg s _I cirT 2 if mergers are taken into account. 



1. Introduction 

Cluster selection by X-ray emission is generally believed to 
produce well-defined samples of massive clusters. This is cer- 
tainly true for relaxed objects near virial equilibrium, but clus- 
ters undergo substantial evolution during the cosmic epoch 
which we can overlook. Numerical simulations demonstrate 
that temperatures and X-ray luminosities of the intracluster gas 
increase by factors of a few for periods which are compara- 
ble to the sound-crossing time while clusters undergo major 
mergers. This may lead to a substantial contamination of X-ray 
flux-selected cluster samples by less massive, but dynamically 
active clusters. This bias needs to be quantified before cosmo- 
logical conclusions based on the cluster population can be con- 
sidered reliable. 

Massive and compact galaxy clusters are also efficient 
strong lenses. This gives rise to the expectation that strong lens- 
ing should be particularly frequent in X-ray selected cluster 
samples, and in fact many X-ray luminous clusters have been 
found to be strong gravitational lenses. 

However, strong clus ter lensing can al so be transiently in- 
creased by factors < 10 dTorri et al. 2004) during major clus- 
ter mergers, on time-scales comparable to the dynamical clus- 
ter time-scale. This may lift relatively low-mass clusters above 
the criticality limit for strong lensing which would otherwise 
be undercritical. Based on a relativ ely small sample of nu- 
meric ally simulated galaxy clusters, iBartelmann & Steinmetz 
(1996) pointed out that X-ray selection is not guaranteed to 
select for the most efficient, strongly-lensing galaxy clusters. 
Early work on the interplay between strong lensing statis- 
tics and observati onal selection effec ts can also be found in 
Wu & Maol d 1996b and lCooravl dl999l) . 
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Both effects of major mergers, the enhancement of their 
X-ray visibility and their strong-lensing efficiency, potentially 
open a huge reservoir of clusters which would remain unob- 
servable in quiescence. The amplitude of this effect must de- 
pend on the frequency of major mergers, and thus on the cos- 
mological model and its parameters. Specifically, merger rates 
at fixed redshift depend on the amount of dark matter and dark 
energy and its cosmic evolution. 

Here, we address the question how X-ray cluster selection 
may affect the strong-lensing efficiency of the selected clusters, 
and what fraction of the optical depth for strong lensing we can 
expect to be produced by galaxy clusters visible above a certain 
X-ray flux limit. In addition to the standard ACDM cosmology, 
we analyse three cosmological models with dynamical dark en- 
ergy which predict different merger histories. Two of them have 
a small, but finite dark-energy density in the early universe and 
are thus expected to host cluster populations which form earlier 
in time. 

We combine tw o semi-analytic methods, one derived by 
Randall et alj 020021) describing the enhancement of X-ray tem- 
peratur es and luminosities during mergers, and another devel- 
oped by lFedeli et al J (120061) for calculating strong-lensing clus- 
ter cross sections. Cluster merger histories are modelled by 
merger trees planted in the extended Press-Schechter formal- 
ism. 

We summarise the model for the enhanced X-ray tempera- 
ture and luminosity in the next section and their translation into 
observable fluxes in Sect. 3. We describe the application of this 
model to simulated cluster populations in Sect. 4, present our 
results in Sect. 5 and conclude in Sect. 6. 
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2. Luminosity and temperature boost 

When galaxy clusters undergo violent dynamical events, such 
as interactions with substantial substructures or mergers with 
galaxy groups and clusters of comparable mass, the intra- 
cluster medium (ICM) is compressed and heated by ram pres- 
sure and shock waves. This results in an overall enhancement 
of the mean gas temperature and the X-ray emissivity due 
to bremsstrahlung. The dynamics of gas and dark matter in 
clusters during major mergers is typically very complicated 
and usually studied bas ed on numerical simulations (see e.g. 
Pool e et al 2006, 2007 for a recent review and applications). 

It will be sufficient for our purposes to model the short- 
term increases in temperature and X-ray luminosity in a sim- 
plified manner which captures their important characteristics in 
a statistically correc t way. Such a simplified model is given by 
Randall et al.1 (120021) . There, the authors derive fitting formulae 
for the time-dependent increase in average temperature relative 
to its unperturbed value of the combined system of main cluster 
and merging body during its interaction. 

They employ A^-body cluster simulations combined with 
adiaba tic hydrodynamics developed and describe d in earlier 
work dRicker & Sarazinll200ll; iRicker etalJ fcoOO). Shocks in 
the ICM are extremely well resolved in these simulations, al- 
lowing temperature and luminosity increases to be studied in 
detail. 



Randall et al.1 d2002l) find that the total time interval At dur- 
ing which the average temperature or bolometric luminosity of 
the ICM of the system are raised above fixed levels T or L is 
given by 



f = V[(r-r c ) 2 -i]( £ 2 + 



(i) 



where F is the ratio between the quantity in question, T or L, 
and its unperturbed value Tq or Lq. The parameter T c is related 
to the maximum value of the boost, as will be explained below. 
The time interval is measured in units of the sound-crossing 
time f sc of the main cluster body by £ = log(Af/f sc ). Assuming 
isothermal gas, the sound-crossing time is 



R v Him v 



(2) 



where the virial radius of the dark matter halo of the main clus- 
ter, R w , is taken as a characteristic dimension, [i is the mean 
molecular weight of the ICM (assumed to be yu = 0.59 through- 
out this work) and m p is the proton mass. 

The merger is characterised by the mass fraction of the sec- 
ondary cluster, 

M 2 



/ = 



M\ + M 2 



which obviously reaches a maximum value of 1 /2 for equal- 
mass mergers. The fit parameters F c and s are then expressed 
as power laws of /, 



and 



s=Cf- 



% C = E [in (Mi +M 2 )-F In (mJ /3 + M^ 3 )] 



(4) 



(5) 



Table 1. Best-fit amplitudes and exponents for the temperature 
(top row) and luminosity (second row) boosts during cluster 
mergers. 



quantity 


A 


B 


C 


D 


E 


F 


temperature 
luminosity 


3.98 
8.28 


0.448 
0.659 


0.96 
0.91 


0.539 
0.316 


3.71 
-0.74 


2.81 
3.29 



The amplitudes and exp onents appearing in the last three equa- 
tions were calibrated bv lRandall et al. (2002) against their sim- 
ulations. They are summarised in TableQ]for the boosts in both 
the temperature and the bolometric X-ray luminosity. 
Solving Eq. (HJ for F, 



r = r c - A i + 



(£-£) 2 



i 



(6) 



we see that it allows a maximum value for the boost of F m = 

r c -i. 

We plot the function T(^) for the temperature and the lu- 
minosity in Figure [1] It shows the temperature and bolometric 
luminosity of the merging system in units of the pre-merger 
values as a function of time in the interval between the begin- 
ning of the boost (r = 1) to the moment of perfect overlap of 
the two clusters (and thus of the maximum boost, F = T m ). The 
mass of the main cluster is set to Mi = 7.5 ■ 1O 14 M , and re- 
sults are shown for four different values for the mass M 2 of the 
merging substructure. 

As intuitively expected, these plots show that the maximum 
temperature and luminosity reached by the system is larger 
when the masses involved in the merger process are similar. In 
that case, duration of the boost is also minimal. Moreover, the 
curves illustrate that the relative increase in bolometric lumi- 
nosity exceeds the one in average temperature, reflecting the 
higher sensitivity of bremsstrahlung emission to the density 
compared to the temperature. 

Finally we emphasise that the fitting formulae for the 
average temperature and bolometric X-ray luminosity en- 
hancements given above are valid only for head-on mergers. 
Generalisations to non-zero impact parameters are given in 
Randall et al.1 (120021) . but we assume head-on mergers through- 
out for simplicity. Note that, for a non-head on merger, the du- 
ration of the boosts in temperature and luminosity is larger, but 
their maximum values are smaller. These are two somewhat 
counter-acting effects, and we expect the total influence to be 
not significant for our pourposes. 



(3 ' 3. Fluxes obtained from individual clusters 



We shall refer mainly to one particular set of observed 
galaxy clusters when comparing to observations, i.e. the 
ROSAT-ESO Flux Limited X-ray c luster sample {Reflex, 



Collins et all 120001 ISchuecker etal] |200U iBohringer etal 



20011) , which was drawn from th e ROSAT All-Sky Survey 
(RASS. ISnowden & SchmitJl 990). In this section, we describe 
the construction of a synthetic cluster sample imitating the pro- 
cedure used for the construction of the Reflex sample. 
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Fig. 1. The enhance ment of temperature (left panel) and luminosity (right panel) of an interacting system according to 
Randall et al.1 (|2002), where the main cluster has mass M\ = 7.5 ■ 1O I4 M , and the substructure's mass takes four differ- 
ent values, as labelled in the plot. Time is measured in units of the sound-crossing time of the main structure and starting from 
the instant to of maximum boost. 



3.1. Ideal flux 

We describe the merger history of individual clusters by means 
of merger trees constructed based on th e extended Press 
Schechter theory (.Press & Schechteri.l974t .Bond et al.|[l991 



Lacey & Coldl993-ISomerville & Kolatdll999h . Thus, the only 
information we have on each individual galaxy cluster is its 
mass and its redshift. We first related these properties to the 
idealised X-ray flux, that is the flux that would be measured 
in the absence of any instrumental issue. Next, we shall add 
background noise, convolution with the point-spread function 
(PSF), and the detector response. 

We start from the virial relation between mass, redshift and 
temperature of the ICM, 



kT = 4.88keV 



M 



10 15 M C 



-Kz) 



2/3 



(7) 



where h(z) is the (reduced) Hubble parameter at the redshift z 
of the cluster, and the nor malisation constant is calibr ated with 
the cluster simulations of Mathiese n & Evrardl (.2001 ). 

Introducing the temperature-mass relation (|7]i into a merger 
tree, we assign temperatures to individual clusters. When a 
cluster is merging with a substructure according to its merger 
tree, we can either ignore the temperature and luminosity boost 
caused by the merger. In this case, only the increasing cluster 
mass will cause the temperature to rise. Or, we can boost the 
temperature according to the description outlined in the previ- 
ous section, depending on the state of the merger process. In 
both cases, we obtain a unique temperature for each cluster at 
each redshift step in its merger tree. 

We note here that no statistical fluctuations are taken into 
account in our assignment of X-ray temperatures and lumi- 
nosities to clusters of a given mass. Thus, once the mass of 
the main cluster body and a merging subclump are fixed, the 



same merger phase will always lead to the same temperature 
increase. However, mergers do introduce statistical fluctuations 
into the temperature-luminosity relations of our simulated clus- 
ters. We shall return to this point further below. 

We also clarify that we consider only binary mergers here. 
Whenever a cluster undergoes a multiple merger, we model 
only the one with the most massive substructure, neglecting 
the others in comparison. Since the simultaneous interaction 
of a galaxy cluster with more than one massive substructure is 
an extremely rare event, we believe that this approximation is 
sufficiently accurate for our pourposes. 

Figure |2] shows the relation between the mass and the tem- 
perature for low-redshift model clusters in our synthetic sam- 
ple, after the temperature boost due to mergers has been ap- 
plied. The normalization is obviously higher (a factor ~ 2) 
than that of the original relation Eq. (0 because cluster merg- 
ers always increase the cluster temperatures. However, the 
slope of the relation is well preserved, and the resulting sam- 
ple fairly reproduces the observ ed HIFLUGCS cluster sam- 
ple dReiprich & Bohringerl.20021) . which mostly contains low- 
redshift clusters. This supports the validity of our model. 

Next, we need to derive the ideal flux from the mass, 
the redshift and the temperature. We do so using the p ubli- 
cally available software package xspec dArnaud 1996), as- 
suming that the ICM can be described by a Raymond-Smith 
(.Raymond & Smith! 1 1977b plasma model. We set the metal 
abund ance to Z _ 0.3Z o (Fukazawa et alj Il998c ISchindlei 
1 19991 see also .Bartelmann & Wh ite 2003). We normalise the 
spectrum by means of the observationally calibrated relation 

/ kT \ 2 ' 331 

^=2.5. 10^^^ ) , (8) 

derived by [Allen & Fabian fll998h. where _-boi jg the bolomet - 
ric X-ray luminosity (see also [Mushotzky & Scharf] (1997); 
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Fig. 2. The mass-temperature relation of our synthetic cluster 
sample. Black, empty circles represent our sample clusters at 
low redshift. The green dashed line is the original M - T rela- 
tion given by Eq. @. The red filled squares with errorbars are 
the observed clusters of the HIFLUGCS sample, and the blue 
solid line is a relation with the same slope as Eq. (0 but with a 
normalization higher of a factor ~ 2. 



Reichartetai1(ll999t) : Hash imoto et a l. (2002) for some discus- 
sion on the redshift e volution of th i s Lt, i - T relation). 

2006t 



Several auth ors dStanek et al 



O'Hara et alj [2006; 



Pratt et al. 2006]) have discussed that the observed scatter in the 
temperature-mass and luminosity-temperature relations might 
not be entirely caused by recent mergers, but rather be sen- 
sitive to the complete merger history of the cluster. However, 
we verified that the scatter in the natural logarithm of mass 
around the best fit mass-luminosity relation that we obtain for 
nearby clusters in our synthetic sample is cn nM - 0.4 and 
thus agrees well with the observed value for the H IFLUGCS 
data dReiprich & Bohringerir2002t IStanek et al1l2006l) . Hence, 
we conclude that our modeling of mergers introduces scatter 
into the luminosity-mass relation compatible with the observed 
scatter, and thus fairly captures the observed statistical fluctua- 
tions. 

Finally, we take account of the interstellar absorption by 
neutral hydrogen in the Milky Way. We do so by combining our 
Raymond-Smith plasma model with the phabs multiplicative 
model component of the xspec software, adopting a constant 
hydrogen column density of n# = 4- 10 20 cm" 2 , appropriat e for 
relatively high galactic latitudes dDickev & Lockmanlll990l) . 



3.2. Instrumental effects 



In con structing the Reflex cluster sample, iBohringer et al 



(1200 ll) used the count rate received for each individual clus- 



ter by the ROSAT PSPC detector in the energy channels cov- 
ering the [0.5, 2.0] keV energy band. To compute synthetic 
count rates for each object in our simulated cluster popula- 
tion, we first need to introduce a model for the distribution of 
the ICM within the clusters. We adopt the isothermal yS-model 
(ICavaliere & Fusco-Femiandl 19761 119 78). assuming B = 2/3 
throughout, following lMohr & Evr ard ( 1999). 
The resulting gas-density profile is 



PicmO) 



PlCM.O 

l+r 2 /r 2 ' 



(9) 



Its core radius r c is related to the X-ray luminosity in the 
[0.5,2.4] keV energy band through 



r c = 0.125Mpc/T 



5 • 10 44 ergs- 



0.2 



(10) 



djones et al.lll998 ; see Vikhlinin et alj[l997 for a discussion on 
the redshift evolution of this relation). 

Since the emissivity of thermal bremsstrahlung is propor- 
tional to the squared gas density, this density profile implies 
the surface brightness profile 



B(d) 



So 



(l+«) 



3/2 ' 



normalised by 



Bo = 



27rt9 2 ' 



(11) 



(12) 



Here, F is the total flux in the respective energy band (i.e. the 
integral of the surface brightness profile over the solid an- 
gle), 6 = r/D A (z) and G c = rJD A (z). Obviously, D A (z) is the 
angular-diameter distance to redshift z. 

We next convolve this surface-brightness profile with 
the instrumental PS F. The shape of the ROSAT -PSPC PSF 
is summarised in i Bartelmann & White (2003) based on 
Hasinger et al. ( 19951 . Its shape depends slightly both on the 
energy channel considered and on the off-axis angle of the 
source. For simplicity, we shall assume on-axis sources and an 
energy channel at 1 keV, approximately at the centre of the en- 
ergy bands considered in the present work. 

Background count rates are provided in form of a map 
on the RASS web pag^H We use a constant median value of 
b = 2.6 • 10~ 4 s _1 arcmin~ 2 , but note that our results are quite 
independent of the background correction. 

Finally, we obtain the count rate produced by the PSF- 
convolved, background-corrected surface-brightness profile. 
We integrate over the complete profile using the f akeit com- 
mand of the xspec software, adopting the PSPC response ma- 
trix in the [0.5,2.0] keV energy band. 



3.3. Nominal flux 

The nominal flux F n for the Reflex cluster sample is defined 
as the flux produced in the [0.1,2.4] keV energy band by a 



1 http : //www . xray . mpe . mpg . de/ cgi -bin/rosat/ 
rosat-survey 
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Fig. 3. The normalized difference between the nominal and the 
ideal flux as a function of the ideal flux itself is shown for the 
synthetic cluster population used in this work. The upper solid 
and lower dash-dotted curves show the mean and the scatter, 
respectively. The nominal flux is always ~ 20% larger than the 
ideal flux, while the scatter in the relation can be neglected only 
for high fluxes. 

Raymond-Smith model plasma set to redshift zero, with a tem- 
perature of 5 ke V, metal abundance of Z = 0.3Z o , absorp- 
tion as given by Dickey & Lockmanl dl990h and a spectrum 
normalised so as to reproduce the observed number counts in 
the energy channels corresponding to the [0.5,2.0] keV band 
dCollins et alJl2000tlBohringer et al.|[200lT) . 

To each cluster and at each redshift step of its merger tree, 
we assign a nominal flux exactly in the same way. The only 
difference between the definitions of our synthetic sample and 
of the Reflex sample is that we normalise the spectrum of the 
plasma model so as to reproduce the count rates computed at 
the end of Subsect. 3.2. 

The Reflex cluster sample is flux-limited, in the sense that 
it contains only galaxy clusters with nominal flux < F n ^ m = 
3.0 • 10" 12 erg s cm 2 . We here adopt the same nominal flux 
as a threshold for synthetic cluster samples. We shall use the 
nominal Reflex flux limit and four additional lower flux limits 
in order to create synthetic samples containing a larger number 
of objects. 

The relation between this nominal flux and the ideal flux 
introduced in Sect. 3.1 (without hydrogen absorption) is shown 
in Fig. [3] where the mean difference between the two fluxes, 
normalised to the ideal flux, and its standard deviation are plot- 
ted as functions of the ideal flux itself. 

According to the definition given at the beginning of this 
section, the ideal flux of a cluster is computed in the energy 
band [0.5,2.0] keV, while the nominal flux is the flux in the 
[0.1,2.4] keV band of a fiducial cluster with fixed physical 



properties (Z = O.3Z , T = 5keV, z — 0) that produces the same 
count rates as the cluster at hand in the [0.5,2.0] keV band. 
Since the nominal flux is computed in a wider and softer band 
than the ideal flux, there is a bias because the nominal flux ex- 
ceeds the ideal flux typically by ~ 20%. The scatter about the 
mean is relatively large for small fluxes, but drops to zero as 
the flux increases because then the effects of PSF convolution 
and background subtraction are smaller. We conclude from this 
plot that the ideal may be used instead of the nominal flux, thus 
saving the time for the computation of hydrogen absorption and 
instrumental effect, but only when the flux is sufficiently large 
(> 10" 12 erg s" 1 crrT 2 ) and accounting for the 20% bias. 

4. Application to the cluster population 

4.1. Cosmological models 

We base our analysis on the s a me co smological models anal- 
ysed in iFedeli & Bartelmannl (12007ft . These are a standard 
ACDM model, a universe with a constant equation-of-state pa- 
rameter Wde = -0.8 for the dark energy, and two early-dark 



energy models (Ferreira & Joyce 1998; Doran et al. 2001a b; 
Caldwell et ail 120031: IWetterichl l2004fK In such models, the 
quintessence energy density tracks that of the dominant com- 
ponent of the cosmic fluid until it begins dominating in the 
late cosmic history. They are characterised by a small but 
non-negligible amount of dark energy at very early times, 
which may influence all kinds of high-redshift processes. 
Bartelmann et al.l (|2006) showed that early dark energy also af- 
fects the non-linear structure formation, giving rise to signif- 
icant differences in the mass function and the merger rates of 
dark-matter halos compared to the concordance ACDM model. 

The matter power spectrum is scale-invariant (n — 1) for 
the models with a constant u>d e , while the spectral indices are 
n = 1.05 in the first and n = 0.99 in the second early-dark 
energy model (hereafter called models EDE1 and EDE2 re- 
spectively). The other cosmological parameters for all mod- 
els studied here are summarised in Tab. [2] They were set 
such as to mat ch the power spectrum o f the CMB temperature 
fluctuations dSpergel et al.l 120031 feOOq). constraints from the 
large-scale structure of the Universe ( Tegmark et al J2 004). and 



observations of type-la supernovae dRiess et alj 2004b . Note 
that this requires different normalisations erg for the dark- 
matter power spectrum in the different cosmological models, 
which also affect the differ e nces i n the mass-assembly his- 
tories. IFedeli & Bartelmannl d2007l) studied the total strong- 
lensing optical depth of the cluster population expected in these 
model cosmologies. Here, we extend this study by taking X-ray 
selection effects into account and exploring the effect of early 
dark energy on the number of clusters in flux-limited samples. 

4.2. Merger trees 

As already stated, we use h ere the merger trees produced in 
Fedeli & Bartelmannl d2007b for a set of N - 500 dark matter 
haloes with masses uni formly distributed between 1 14 M B ft 



and 2.5 • 10 15 M o /z _1 . In Fed eli & Bartelma nn (2007), we com- 



puted for each cluster the lensing cross section for gravitational 
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Table 2. Parameters of the four cosmological models studied 
here. 





EDE1 


EDE2 


w dc = -0.8 


ACDM 




0.82 


0.78 


0.80 


0.84 


h 


0.67 


0.62 


0.65 


0.65 


"m,0 


0.33 


0.36 


0.30 


0.30 




0.67 


0.64 


0.70 


0.70 



arcs with length-to-width ratio larger t han d — 7.5, usin g the 
fast, semi-analytic method developed bv lFedeli et al.l2 006. The 
objects in the sample were evolved up to a source redshift ran- 
domly drawn from an observed distr i bution of faint blue back- 
groun d galaxies dSmail et al.l 1 19951: iBartelmann & Schneider! 
20011) . Lensing efficiencies were calculated both ignoring and 
accounting for mergers with substructures which transiently 
enhance the cross sections. 

We apply the same scheme here, calculating the nominal 
X-ray flux for each sample cluster at each redshift step. We 
do that both ignoring and accounting for the effect of mergers 
which transiently enhance the intrinsic luminosity and temper- 
ature (and thus also the nominal flux) of the clusters. 

Finally, we can combine this information with the calcu- 
lations of the lensing efficiencies, evaluating the effect of flux 
selection on the statistics of gravitational arcs in galaxy clus- 
ters. 

4.3. Cluster number counts 

It is interesting note in passing how the previous calculations 
predict the total number of clusters in a flux-limited sample to 
change as a function of the limiting flux in the different cosmo- 
logical models, and what the quantitative effect of mergers is in 
this respect. 

We emphasise here that our algorithm for producing 
synthetic cluster samples is not ideally adapted to cluster- 
abundance studies because our cluster sample includes only 
relatively high-mass haloes, and thus the low-mass end of the 
distribution is not well sampled. Moreover, we did not take into 
account the likely steepening of the temperature-luminosity re- 
lation for low-mass clusters or galaxy groups. While this has 
no effect on samples of X-ray luminous, hot or strongly lensing 
clusters, it is likely that the overall number of structures in the 
different flux-limited samples is overestimated. This is because 
we tend to assign to objects with very low mass a temperature 
higher than expected in presence of steepening. 

In Figs. [Hand [5] we show the total number of galaxy clus- 
ters predicted to be observed in the four cosmological models 
used here as a function of the limiting nominal flux. We plot 
results obtained by accounting for and ignoring the effects of 
cluster mergers, and indicate the total number of galaxy clus- 
ters observed in the Reflex cluster sample, extrapolated to the 
whole sky. 

Several interesting pieces of information can be read off 
these figures. First of all, cluster mergers increase the total 
number of visible objects by factors between 2 and 3. This fac- 





L °g( F n,iim[ er g s_1 cm" 2 ]) 

Fig. 4. The total number of galaxy clusters observable on the 
whole sky, given as a function of the limiting nominal flux for 
the four different cosmological models considered here. Black 
solid lines are obtained by ignoring the transient boost due to 
cluster mergers, red lines are obtained taking it into account. 
The cyan horizontal line gives the number of galaxy clusters 
obtained from the Reflex cluster sample extrapolated to the 
whole sky (F n ,iim = 3 ■ 10~ 12 erg s _1 cm -2 ). 



tor tends to decrease towards lower flux limits because of two 
effects. First, at low flux limits, the total number of clusters ob- 
servable without mergers is larger, thus the fractional increase 
due to cluster interactions tends to be smaller. Second, at low 
flux limits, we include low-mass objects into the sample whose 
merger frequency is lower. We also see that, according to this 
analysis, only the models with constant equation-of-state pa- 
rameter are in agreement with the Reflex observations, while 
early-dark energy models overpredict the cluster abundance by 
a factor of « 2. 

In Fig. [6] we fix the nominal flux limit to that of the Reflex 
sample, F n n m = 3 • 10~ 12 erg s crrT 2 and show a histogram of 
the redshift distribution of observed clusters. Again, we show 
results with and without the enhancements by cluster merg- 
ers. For all models, the number of clusters drops to zero above 
z ~ 0.3 with mergers, and already above z ~ 0.15 without 
mergers. The absence of substantial differences between differ- 
ent cosmologies is due to the fact that at low redshift the differ- 
ence between the structure formation in presence or absence of 
early-dark energy tends to d isappear (see also the discussion in 



Fedeli & Bartelmann 2007). It is interesting to note that the ob- 



servational results from the Reflex sample dCollins et al . 2000) 
are qualitatively very well reproduced only accounting for clus- 
ter mergers. Ignoring the effect of interactions in our models 
leads to an underestimate of objects in the high-redshift tail. 

In the context of Fig. [6j we also note that while the quali- 
tative trend and the peak position of the observed Reflex distri- 
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Fig. 5. As Fig. |U but combining all curves in the same plot. 
Black and red curves are for the EDE1 and EDE2 models, 
respectively. The blue curve is for the model with constant 
equation-of-state parameter Wd e = -0.8, and the green line is 
for the ACDM model. For each model, the solid and dashed 
curves are obtained ignoring mergers and taking them into ac- 
count, respectively. The horizontal line shows the number of 
clusters observed in the Reflex cluster sample. 

bution are reproduced, the normalisation is generally too high, 
for the reasons discussed. 



5. Results 

We finally return to the main purpose of this work, that is probe 
how cluster selection by their X-ray flux may influence the 
optical depth of the sample for the production of pronounced 
gravitational arcs when the effect of cluster mergers are taken 
into account. In other words, we analyse how the total num- 
ber of arcs that we can expect to observe in an X-ray selected 
galaxy-cluster sample depends on the X-ray flux limit of the 
sam ple itself. While we have co nsidered only arcs with d > 7.5 
in 



Fedeli & Bartelmannl (12007b . we extend the analysis here to 



arcs with d > 10 (so-called giant arcs). 

Figure [7] shows contour lines in the mass-redshift plane 
for the nominal X-ray flux of clusters from our synthetic sam- 
ple, and for the cross section for arcs with length-to- width ra- 
tio larger than d = 7.5 for the ACDM cosmological model. 
Contours in the left and right panels were obtained ignoring 
mergers and taking them into account, respectively. We over- 
plot the flux limit for the Reflex cluster sample. The edge in 
both the X-ray flux and lensing-efficiency contour lines going 
from the upper left to the lower right corner illustrates the lack 
of high-mass clusters at high redshift. These figures clearly 
show the effect of mergers on the lensing efficiency and the 
average X-ray flux from clusters. Including mergers, the con- 




0.4 

Redshift 

Fig. 6. Histograms for the cluster redshift distributions ex- 
pected in the cosmological models studied here, as labelled in 
the boxes. The limiting flux is the same as for the Reflex cluster 
sample, F n ^ m = 3 • 10~ 12 erg s _1 cm 4 . As before, the black 
lower curve does not account for cluster mergers, while the red 
upper curves do. 



tour lines are much more irregular and extend towards lower 
masses and higher redshifts, both for the nominal X-ray flux 
and for the cross section. The geometric suppression of the 
lensing efficiency at very low redshift is also evident since the 
black contours in the lower panels never reach z = 0, and turn 
further away from z = for lower cluster masses. The lens- 
ing efficiency expected to be observed in a ACDM model in 
a Reflex-like cluster sample is thus contributed only by those 
clusters falling between the cyan curve and the lower contour 
lines in the lower panels. 

We emphasise that Fig.[7jonly shows the properties of indi- 
vidual objects in our synthetic sample. In order to find sample 
properties, this information must be convolved with the cluster 
mass function. This means that even the small wiggles appear- 
ing in the heavy contour in Fig. [7] when mergers are included 
will have a substantial effect on lensing statistics and number 
counts. We prefer not to weight with the mass function in Fig. [7] 
to illustrate exclusively the effect of the flux cut. To have an 
idea of the consequence that small shifts in the M — z plane 
can have on the quantitative results, the number of structures in 
the narrow mass interval 14.5 < logM < 14.6 doubles when 
computed for z < 0.10 compared to z < 0.13. 

Figure [8] shows the average optical depth for gravitational 
arcs with length to width ratio d > 7.5 and d > 10 pre- 
dicted to be observed in a flux-limited X-ray cluster sample as 
a function of the limiting flux. Results obtained both ignoring 
and taking account of clus ter mergers are shown. As noted in 
Fedeli & Bartelmannl ( 2007 ) before, mergers increase the opti- 
cal depth by a factor between 2 and 3. The present figure shows 
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Fig. 7. Contour lines for the nominal flux (upper panels) and for the lensing cross section for gravitational arcs with d > 7.5 
(lower panels) in the ACDM model. The thick cyan contours correspond to the limiting nominal flux of the Reflex cluster sample. 
The upper heavy black contour line in the upper panels corresponds to a flux of 1CT 14 erg s _1 cm 2 , the others are spaced by one 
order of magnitude. Analogously, the lowest black contour in the lower panels corresponds to a cross section of 10~ 4 Mpc 2 /z~ 2 , 
and the others are spaced of one order of magnitude. Merger processes are taken into account in the right panels, while they are 
ignored in the left ones. 



that this remains true for all flux limits considered here even 
when X-ray selection effects are taken into account. Moreover, 
we note that the slope of the f - F n \im relation tends to increase 
(decrease in absolute value) towards low limiting fluxes. This 
is due to the fact that the lensing efficiency drops above z ~ 0.3, 
and approaches zero towards sufficiently high redshifts. Thus, 
if the flux limit is low enough, the sample contains all the arcs 
that are produced and that would be observed without selection 
effects. 

Increasing the length-to-width threshold from d = 7.5 to 
d = 10 changes the absolute value of the average optical depth, 



but not its qualitative behaviour as a function of the limiting 
flux. 

A central result of our study is the ratio between the opti- 
cal depth for large gravitational arcs in an X-ray selected, flux- 
limited cluster sample compared to the total, idealised optical 
depth. Since the average optical depth is related to the num- 
ber of arcs by a constant factor (determined by the total num- 
ber density of background sources), this corresponds to the ra- 
tio between the number of arcs expected to be seen in a flux- 
limited sample of X-ray clusters, and the total number of arcs 
that would be observable in absence of any X-ray selection ef- 
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Fig. 8. The observed average optical depth for arcs with length to width ratio larger than d = 7.5 (left panel) and d — 10 
(right panel) as a function of the nominal X-ray flux limit. Results are shown for the four different cosmologies considered here. 
Black solid lines represents results obtained by ignoring mergers processes, while broken red lines are obtained taking them into 
account. 



feet. This ratio is shown in the right and left panels of Fig. [9] 
accounting for and ignoring the transient merger boosts in tem- 
perature, X-ray flux and lensing cross section, respectively. 

We note that, when we include mergers, the fraction of the 
average optical depth in an X-ray selected cluster sample in- 
creases with increasing limiting flux. We attribute this to the 
facts that (i) cluster mergers tend to enhance the lensing effi- 
ciency, which affects both the X-ray selected and the total av- 
erage optical depth, and (ii) cluster mergers also enhance the 
clusters' temperature and the X-ray flux, which only affect the 
X-ray selected optical depth. We also note that there is a slight 
tendency for models with lower power-spectrum normalisation 
o"8 to have a larger fractional importance for the observed av- 
erage optical depth. For instance, the ratio between the optical 
depth of X-ray selected clusters to the total optical depth is 
systematically slightly larger for the model EDE2 (<x 8 = 0.78) 
than for the model EDE1 (erg = 0.82). However, this effect is 
very small. 

6. Summary and conclusions 

We have studied here the influence of selection effects on the 
total observed number of gravitational arcs in X-ray selected 
galaxy-cluster samples, taking cluster mergers into account. 

To perform our study, we considered the assembly history 
of a synthetic sample of galaxy clusters. We linked the virial 
mass of the dark-matter cluster halo to the temperature of the 
ICM using the virial relation (0. We then used t he analytic fit- 
ting formulae provided by Rand all et al. I (120021) for the boost 
in temperature caused by merger processes between clusters 
and substructures during the formation. After wards, we used 
the publicly available software package xspec dArnaudl 19961) 
to convert the (boosted and unboosted) temperature of the ICM 



into the ideal flux (in front of the instrument) produced by each 
individual object accounting for redshift, metal emission lines 
and interstellar absorption. Then, using the response matrix of 
the PSPC detector on-board the Rosat satellite, we transformed 
the ideal flux in a photon count rate, also including observa- 
tional effects such as realistic background count rates and PSF 
convolution. We finally turned the count rates for each ob- 
ject into a nominal flux, as defined in the construction of the 
ROSAT-ESO Flux Limited X-ray cluster sample Reflex. 

We repeated this procedure to the four different cosmo- 
logical models previously used in Fedeli & Bartelmannl ( 2007 ) 
for studying the lensing properties: A ACDM model, a dark- 
energy universe with constant equation-of-state parameter w = 
-0.8 for the dark-energy component, and two early-dark en- 
ergy models, whose relevance in terms of linear and nonlinear 
structure formation has recently been pointed out. 

As an intermediate, qualitative result, we analysed the to- 
tal number of clusters expected to be visible in X-ray selected 
cluster samples as a function of the nominal flux limit. We ob- 
tain a significant difference between cosmologies with early- 
dark energy and the model with constant equation-of-state pa- 
rameter, and also a significant difference due to the introduc- 
tion of the effect of cluster mergers. In particular, we find that 
the qualitative redshift distribution of clusters observed in the 
Reflex sample can only be reproduced accounting for merger 
boosts in temperature and luminosity. Moreover, early dark- 
energy models seem to overpredict the total number of ob- 
served objects, although our results are not precise at the low- 
mass end of the distribution, which is irrelevant for strong 
lensing. Thus, our absolute numbers are likely to be an over- 
estimate. This may hint at a potentially very interesting test 
for early-dark energy in particular, and for the dynamics of 
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Fig. 9. The fraction of the optical depth for the formation of large gravitational arcs, contributed by X-ray selected clusters in a 
flux-limited sample, compared to the total optical depth. Different line styles represent different cosmological models, as labelled 
in the plot. The effects of cluster mergers are taken into account in the right panels and ignored in the left panels. The upper and 
lower panels show results for arcs with length-to-width ratios exceeding d = 7.5 and d > 10, respectively. 



quintessence models in general, and certainly warrants further 
investigation in the future. 

Finally, our main results are determined by the combination 
of the nominal X-ray flux selection with the observed strong- 
lensing statistics. We find that cluster mergers enhance the av- 
erage observed optical depth by factors between 2 and 3 for 
all limiting fluxes considere d here. We also confirm the result 
Fedeli & Bartelmann (2007) obtained ignoring any X-ray se- 
lection effects, that the different structure-formation history in 
early-dark energy models causes the lensing efficiency to in- 
crease by a factor of * 3 compared to models with constant 
equation-of-state parameter for the dark energy. This remains 
true for all limiting fluxes, and we see that, for instance, the 



same lensing efficiency reached in a ACDM model with the 
help of cluster mergers is reached ignoring mergers in mod- 
els with early-dark energy, because of the higher cluster den- 
sity at moderate and high redshifts. Moreover, the slope of the 
f a -F a a m relation tends to flatten towards lower limiting fluxes, 
indicating that we are approaching the total average optical 
depth. 

In an earlier study, iMeneghetti et al.1 J2005) addressed the 
difference between the total strong-lensing optical depths of 
cluster populations in cosmological models with different types 
of dark energy, finding that dynamical dark energy tends to in- 
crease the optical depth for giant arcs compared to the stan- 
dard ACDM. While this seems to be a generic prediction for 
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dark-energy models, they used numerical simulations ignoring 
X-ray selection effects, while we semi-analytically model the 
cluster population, its lensing properties and selection issues. 
We extend our method to models with an early dark-energy 
contribution. 

We also find that the ratio of the flux-limited to the total 
average optical depth is larger when we consider the effect of 
cluster mergers than when we ignore it. This is due to the fact 
that the increment in the lensing efficiency affects both the flux 
limited and the ideal optical depth, while the boost in tempera- 
ture and luminosity affects only the former. 

We carried out these calculations for gravitational arcs with 
length-to-width ratios d > 7.5 and d > 10. We find an (ex- 
pected) difference in the absolute value of the average optical 
depth, while the trend with the limiting X-ray flux is qualita- 
tively unchanged. 

Predicting the number of long and thin gravitational arcs 
to be observed in X-ray selected cluster samples and in dif- 
ferent cosmological models will be very useful in the near fu- 
ture. Forthcoming strong lensing surveys ( Cabanac et alj |2007) 
and the development of automatic detection algorithms for 
strong lensing features dLenzen et alj2004 ; Horesh et al.l l2005 



ISeidel & Bartelmann! l2006h will allow to place further con- 
straints on the dynamics of structure formation in a universe 
dominated by dark energy. 
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